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As cell proliferation is limited due to the loss of telomere repeats in DNA 
of normal somatic cells during division, telomere attrition can possibly play 
an important role in determining the maximum life span of organisms as well 
as contribute to the process of biological ageing. With computer simulations 
of cell culture development in organisms, which consist of tissues of normal 
somatic cells with finite growth, we otain an increase of life span and life ex- 
pectancy for longer telomeric DNA in the zygote. By additionally considering 
a two-mutation model for carcinogenesis and indefinite proliferation by the 
activation of telomerase, we demonstrate that the risk of dying due to cancer 
can outweigh the positive effect of longer telomeres on the longevity. 
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^ ■ 1 Introduction 
o 

Telomeres are tandem repeated noncoding sequences of nucleotides at both ends of 
the DNA in eukaryotic chromosomes stabilizing the chromosome ends and prevent- 
ing them from end-to-end fusion or degradation [ij. Polymerase cannot completely 
replicate the 3' end of linear DNA, so telomeres are shortened at each DNA replica- 
tion [2]. This end replication problem leads to a finite replicative capacity for normal 
somatic cells They can only divide up to a certain threshold, the Hayflick limit 
(HIH]. The enzyme telomerase, repressed in most normal somatic cells, synthesizes 
and elongates telomere repeat sequences at the end of DNA strands so that certain 
cells like germline cells are immortal and indefinite in growth |H1IZ|- 

Most forms of cancer follow from the accumulation of somatic mutations |E1 03 . 
Cancer-derived cell lines and 85-90% of primary human cancers are able to synthesize 
high levels of telomerase and thus are able to prevent further shortening of their 
telomeres and proliferate indefinitely [10J. But if cells are premalignant or already 
cancerous and telomerase is not yet activated, the proliferative capacity of these 
cells and therefore the accumulation of mutations is determined by the remaining 
telomere length So the frequency of malignant cancer should be higher for 

longer telomeres in normal somatic cells. 

Recently published data show that longer telomeric DNA increased the life span 
of nematode worms \V2\ . Experiments are running ^31 whether there is a positive 
effect on the longevity also for organisms with renewing tissues if the telomere length 
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in fetal cells is increased [13j. As the probability for the incidence of cancer is 
correlated to the replicative potential of the mutated cells [Hj, one can ask the 
following question: Is an extension of life span possible if telomeres in embryonic 
cells are elongated and cancer is considered? An answer to this question could be 
given by computer simulations as the presented model focuses on organismal ageing 
due to the loss of telomeres in DNA and neglects other effects which lead to a 
decreasing survival probability with age. 

As shortening of telomeres is one of the supposed mechanisms of ageing on cellular 
level, most stochastical and analytical studies investigate this relationship [H3 HE! 
[13 [TBI EH- A theoretical model which directly relates telomere attrition to human 
ageing was first suggested by Aviv et al. [20j. Here we present a different telomere 
dynamics model providing requirements to study the effects of (a) different mean 
telomere lengths in constituting cells, as well as of (b) somatic mutations leading 
to cancer progression and indefinite proliferation due to telomerase activity on the 
ageing of organisms. 

2 Basic model of biological ageing due to telomere 
attrition 

Every organism is developed from a single progenitor cell, the zygote (figure H}. 
The initial telomere length of zygote cells is assumed to be normally distributed with 
mean \i z and standard deviation a z [21]. Telomere repeats lost per division (TRLPD) 
are randomly chosen at each division of every cell from a normal distribution with 
mean Utrlpd arid standard deviation (Ttrlpd- A dividing cell produces a clone who 
inherits the replicative capacity of the progenitor cell at this age. Cells can divide 
until nearly all their original telomeres are lost. 

For every organism the dynamics of the model is as follows: Divisions of the 
zygote and the stem cells derived from it occur 6 times in the early embryo. Each of 
these cells is the progenitor of one tissue. This is followed by a period of population 
doublings where all cells divide once in every timestep until 2 7 cells are present in 
each tissue. In the following maturation stage, cells are chosen randomly for division 
until each tissue reaches the adult size of 10 4 cells. It takes about 26 timesteps until 
an organism is mature. 

Ageing starts now. In every timestep first cells die with 10% probability due to 
events like necrosis or apoptosis. 10% of the cells of the corresponding tissue are 
then randomly chosen for division to fill this gap. The replacement does not have to 
be complete as the chosen cell could probably not divide anymore due to telomere 
attrition. After some time the tissue will start shrinking. The random choice of 
dying and dividing cells in differentiated tissues is in accordance with nature as 
for example in epithelium the choice of cells to be exported from the basal layer is 
random [2ll[2H]. The organism dies, if its total cell population size shrinks to 50% of 
the mature size. The results presented in the following do not depend qualitatively 
on the choice of this threshold. 





Ageing 



Tissue - 10% of all cells die in every timestep due to wear and tear, necrosis, etc. 

^10% new cells fill this gap. Cells are randomly chosen for division, 
but some cells cannot divide anymore due to telomere attrition. 

»- Tissues shrink 



The organisms dies if its size shrinks to 50% of its mature size 



Figure 1: Dynamics of cell proliferation in one organism in the basic model 



3 Results without cancer 

In the implementation of this model, linear congruential generators are used to 
produce the required random numbers and normally distributed variables are gener- 
ated with the Box-Muller algorithm [221 E3| • Resulting age distributions for different 
mean telomere lengths in the zygote cells are shown in (figure |2|). The shape of these 
distributions is analogous to empirical data of many human and animal populations. 
We obtain a positive effect on the longevity of the organisms if the mean telomere 
length in the precursor cell is increased. The chosen mean doubling potentials of 
the zygote cells are 30, 40 and 50 with the choice of \l z = 1500, 2000, 2500 and 
I^trlpd = 50. The number of mitotic divisions observed in human fibroblasts is 
higher [22], but the choice of this parameter is reasonable as the number of con- 
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Figure 2: Age distribution of 32000 organisms with telomere lengths of fj, z 
1500(+), fj, z = 2000(x) and fi z = 2500(*); a z = 100, ^trlpd = 50, a TRLPD = 10 



sidered cells per mature organism (640000) in this model is also much lower than 
in human organisms where the total number of cells is of the order of 10 13 [24J. 
Non-dividing cells are not included in the model as we focus on ageing due to the 
progressive shrinking of tissues driven by telomere shortening. Our results will be 
given below as part of Fig. 5. 



4 Introducing carcinogenesis and telomerase 

Clonal cancer is now introduced in the model. In accordance to the model of Mool- 
gavkar et al., one of our assumptions is that malignant tumors arise from indepen- 
dently mutated progenitor cells (2jj- For most forms of carcinoma, transformation 
of a susceptible stem cell into a cancer cell is suggested to be a multistage process 
of successive mutations with a relatively low probability for the sequential stages 
[HIEHl- Two independent and irreversible hereditary mutation stages are considered 
here, which can occur at every level of development of the organism during cell 
division. 

The first premalignant stage to be considered is a promotion stage: A dividing 
cell can mutate with small probability p mut [24J. All descendant cells inherit this 
mutation. This mutation leads to a partial escape from homeostatic control of 
growth by the local cellular environment (TTJ EH| - Cells on the promotion stage have 
a selective advantage over unaffected cells [HOI- In our model they are chosen first 
for division during maturation and for filling up the gap in the ageing period. 

The subsequent transition can occur again with probability p mu t during division. 
If a cell reaches this second stage of mutation it is a progenitor of a carcinoma. An 
explosive clonal expansion to a fully malignant compartment happens [213 • This cell 



and the clonal progeny doubles in the current timestep until it is no more possible 
due to telomere attrition. This expansion leads only to an increase of the malignant 
cell population size. As a certain fraction of cells is killed per unit time and clonal 
expansions only occur with a very small probability, the tumor environments may 
not continue growing, eventually shrink, or even die j^I]. We assume that it is 
necessary for advanced cancer progression and therefore for the development of a 
deadly tumor that fully mutated cells are able to activate telomerase [32]. 

In our model, telomerase activation is possible at every age of the organism in 
normal and mutated cells during division with a very low probability pteio- The 
irreversible loss of replicative potential is stopped in these cells. As the contribution 
of telomerase to tumorigenicity is not yet completely understood [HHl EE], we 
assume that death of an organism due to cancer occurs if telomerase is reactived in 
at least one fully mutated cell |3H!- We treat the time interval between the occurence 
of the deadly tumor and death as constant, so we set this interval to zero. 

5 Effects of different telomere lengths considering 
cancer 

Figure El shows simulation results for \i z = 1500 and fi x = 2500. As we consider a 
35000 i , , , , ■ , , , , 1 
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Figure 3: Age distribution of 32000 organisms with telomere lengths of \i z = 1500 
with (*) and without cancer (+) and \i z = 2500 with (□) and without cancer 
(x); a z = 100, Utrlpd = 50, (Jtrlpd = 10. Cancer mutations are possible with 
Pmut = 5 * 10 ^. Telomerase can be activated with pteio = 

10" 5 . 

lower complexity by chosing a lower number of tissues and cells per organism, we 
assume higher mutation rates for the incidence of cancer than observed in nature 
[H3EH1- The age distribution for shorter initial telomere lengths considering cancer 



is shifted to the left but still very old organisms exist. For longer telomeres the age 
distribution is again shifted to the left but even behind the distributions for shorter 
telomeres with and without considering carcinogenesis. 
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Figure 4: Mean age of death for different inital zygote telomere lengths jj z (32000 
organisms) without cancer (x) and introducing carcinogenesis (+). 

Thus without considering cancer, organisms with longer zygote telomeres live 
longer, as the life expectancy of the organisms increases linear for longer telomeres. 
But if cancer is considered this effect is reversed for longer initial mean telomere 
lengths (fig. HI). 

The force of mortality P9J resulting from this model is shown for fi z = 1500 with 
and without considering cancer in comparison to empirical human mortality data 
(figure [HI). It agrees to some extent with human mortality functions provided cancer 
is incorporated into the model and decelerates at advanced ages, as claimed for 
human and animal populations |40l 141) . The hump in the curve at younger ages, 
occurring also for other parameter sets, fits to data of many human mortality tables. 



6 Conclusion 

The expected simulation result of the basic model without cancer is an increase of 
life span of most organisms with longer initial telomeres. After introducing somatic 
mutations promoting cancer and telomerase activation in this model, the survival 
probability is lower for each considered initial telomere length in certain time inter- 
vals in adult ages. 

But even low probabilities for the two mutation stages and for the activation of 
telomerase lead to a strong reduction of life span for longer telomeres. So the impli- 
cation of two-stage carcinogenesis for the incidence of cancer in this simple model 
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Figure 5: Mortality function for /i 2 = 1500 with (x) and without cancer (+), 32000 
organisms considered, a z = 100, ^trlpd = 50, utrlpd = 10; German men (*) from 
www.destatis.de (June 2004) (Sterbetafel 2000/2002) 



of cell proliferation in organisms is that life expectancy and life span of complex 
organisms cannot be increased by artificially elongating telomeres in primary cells, 
for example during a cloning procedure. 

Further improvements, extentions and applications of this model are possible. 
With respect to the role of telomeres and telomerase in carcinogenesis, maybe this 
computational approach can contribute to the development of a comprehensive the- 
oretical model in oncology uniting mutagenesis and cell proliferation |4~2] . 
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